Relaxation of a high-energy quasiparticle in a one-dimensional Bose gas 



O 

(N 

as: 
< 



C/2 

c3 



c 

o 

> 

o 

m 

O 

(N 
O 
O 



X 



Shina Tan, 1 Michael Pustilnik, 2 and Leonid I. Glazman 1 

1 Department of Physics, Yale University, New Haven, CT 06520, USA 
2 School of Physics, Georgia Institute of Technology, Atlanta, GA 30332, USA 

We evaluate the relaxation rate of high-energy quasiparticles in a weakly-interacting one- 
dimensional Bose gas. Unlike in higher dimensions, the rate is a nonmonotonic function of tem- 
perature, with a maximum at the crossover to the state of suppressed density fluctuations. At the 
maximum, the relaxation rate may significantly exceed its zero-temperature value. We also find the 
dependence of the differential inelastic scattering rate on the transferred energy. This rate yields 
information about temperature dependence of local pair correlations. 
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Recent experiments with ultracold atomic gases 0, Q 
renewed the interest in fundamental properties of the el- 
ementary excitations in interacting Bose systems. 

A three-dimensional (3D) Bose gas undergoes the 
Bose-Einstein condensation (BEC) phase transition at 
a sufficiently low temperature [3J. The transition af- 
fects dramatically the spectrum of elementary excitations 
(quasiparticles) of the system. In the Bose-condensed 
phase, the quasiparticles obey the Bogoliubov dispersion 
relation e q = sq^/l + (q/2ms) 2 which interpolates be- 
tween a phonon-like linear spectrum at small momenta 
(here s is sound velocity and m is each boson's mass) and 
a free-particle-like spectrum at large momenta [3j. 

The BEC transition affects strongly the lifetime of low- 
energy quasiparticles. The relaxation rate T q of quasipar- 
ticles in the phonon part of the spectrum is very sensitive 
to both their momenta q [4] and temperature T 043; 
Tq oc max{(e 9 ) 5 , e 9 T 4 }. However, the relaxation rate of 
high-energy quasiparticles is dominated by collisions with 
large momentum transfer, does not depend on either q or 
T, and thus is not sensitive to BEC transition Q. Some 
of these long-standing predictions have been recently ver- 
ified experimentally, see [l], 0| for a review. 

Unlike its 3D counterpart, the one-dimensional (ID) 
interacting Bose gas turns at low temperatures to a qua- 
sicondensate in which the long-range order is destroyed 
by quantum fluctuations 0, and the BEC transition 
turns to a crossover. Yet, despite this difference, the spec- 
trum of elementary excitations in ID is still described 
very well by the Bogoliubov dispersion relation . 

However, the quasiparticle lifetime in ID is very dif- 
ferent from that in higher dimensions and is not as 
well understood. The reason is that, due to the con- 
straints imposed by the energy and momentum conser- 
vation, two-particle collisions do not lead to a relaxation 
in ID. At the same time, realizations of ID Bose systems 
with cold atoms confined in tight atomic waveguides [l| 
are described rather well @ by a model of bosons with 
zero-range repulsive interaction (the Lieb-Liniger model) , 
which is integrable [1, U3] ■ In this model, the redistri- 
bution of the momenta between particles in a collision, 
and, therefore, relaxation, are absent 10]. Such apparent 



lack of relaxation was recently demonstrated experimen- 
tally ll| (see the discussion below). 

The leading corrections to the Lieb-Liniger model have 
the form of a local three-particle interaction term [12, EH > 
which breaks the integrability and brings about the 
quasiparticle relaxation. In this Letter, we study relax- 
ation of a particle with a large momentum. This problem 
was considered recently in |13j | , where the inelastic relax- 
ation rate due to three-particle collisions was evaluated 
in the approximation that neglects two-body repulsion. 
The results of Ref. [l3| suggest that, very much like in 
3D, the relaxation rate of high-energy quasiparticles is in- 
dependent of momentum and temperature. However, in 
the present Letter, we demonstrate that, in a dramatic 
departure from the behavior in higher dimensions, the 
relaxation rate in ID depends strongly on temperature 
even at large momenta. It has a pronounced peak at the 
crossover to the quasicondensate state. 

We evaluate the differential and the total relaxation 
rates. Both can be inferred from observations of colliding 



clouds of cold atoms 11, 1 



To describe the relaxation in a weakly-interacting ID 
Bose gas, we consider the simplest Hamiltonian 



H = H a + V, 



(1) 



where 



H =Jdx^(x)(-^-^pjij(x)+^Jdx:p 2 (x); (2) 

describes ID bosons with a repulsive contact interaction 
(hereinafter we set ks = h = 1), and 



V = -^j d x: P \x): 



(3) 



represents the leading integrability-breaking perturba- 



tion 



12j, ll3j. In Eqs. © and ©, p(x) = ^{x)ip{x) 



is the local density operator and the colons denote the 
normal ordering. The strength of the interaction [repre- 
sented by the second term in Eq. ((2J] is characterized [|| 
by the dimensionless parameter 7 = mc/n, where n is 



2 



the ID concentration. A finite three-particle scattering 
amplitude appears already in the first order in a <C 1. 

In this Letter, we study relaxation of a boson with 
momentum q (we assume that q > 0) and kinetic energy 
£ 9 = q 2 /2m, which is large compared to both temperature 
T and a typical interaction energy per particle ui s , 



fq > m&x{T,w s }, 



ms 2 /2. 



(4) 



In the limit of a weak interaction 7 < 1, which we con- 
sider from now on, the sound velocity s in Eq. (j4]) is given 
by [H s = (71/771)^/7. The condition (|4|) ensures that the 
particle is added to an almost empty single-particle state: 

/ 9 = (^)<i ESQ. 

In the lowest (second) order in a the differential rate 
of inelastic scattering is given by 



a q (u) 



27T777/ 



9/3 



dp5{u)-£ q + £, q - p )Q{p,w), (5) 



where Q(p,ui) — jdxdte 1 ^ 1 lpx Q(x,t) is the Fourier 



iwt— ipx 

transform of the correlation function 



g(x,t) = (:p 2 (x,t)::p 2 (0,0):), 



(6) 



which should be evaluated for the Lieb-Liniger model 
Eq. ([2]). In writing Eq. ([5]) we took into account the 
kinematic constraint p < q/3 on the momentum transfer 
in the course of three-particle scattering. The constraint 
translates into a restriction on the transferred energy: 
(J q (to) vanishes for ui > 5£ 9 /9. In terms of a q (u>), the 
total relaxation rate is given by 



r n 



did <J q {bS). 



(7) 



The differential rate (JSJ at large energy transfer ui is 
determined by the behavior of Q(x, t) at t —> 0, 

g(x,t) = 2n 2 g 2 * 2 , *( M ) = (-I^) 1/2 e »- 2 /2*. (g) 

Here r?gi = (: / o 2 (0,0):) is the probability of finding two 
bosons at point x — at time t = 0, and ^(x,t) is 
the solution of the single-particle Schrodinger equation 
with the initial condition ^(x,0) = S(x). Interactions 
do not affect the time evolution in Eq. © as long as 
\t\ <C min{l/o; s , 1/T}. Instead, the dependence on tem- 
perature and on the interaction strength enters Eq. © 
via the normalized local pair correlation g 2 . For the Lieb- 
Liniger model this quantity can be evaluated exactly [13] ■ 
For a weak interaction gi increases monotonic ally with T 
from 52 = 1 at T < T s to g 2 = 2 at T > T s (Tj, where 
we introduced a characteristic temperature scale 

T s = vW?o = ras ; (9) 

here To = 2n 2 /m is the quantum degeneracy tempera- 
ture. (Note that ui s <C T s <C To for a weak interaction.) 



Substitution of Eq. (|5J) into Eq. ([5]) yields the differen- 
tial rate at large positive energy transfer |l8[ 



<T q (u) 



a 2 T g 2 



1 + y/l- oj/£ q 



(l-u/^)(l + 3y/l-u/^) 



1/2 



(10) 

Eq. (ITU1) is applicable at max{cj s , T}«w< 5£ 9 /9. Away 
from the upper end of this interval, at lu <C £ 9 , Eq. (fTU)) 
reduces to 



a a [UJ 



a 2 T Q g 2 



2ir^/W q 



(11) 



To further analyze <J q {ui) at \u>\ -C £ g , we note that 
in this range of ui the momenta p contributing to the 
integral in Eq. ([5]) are small, \p\ -C 9, and it simplifies to 



<j q (w) 



2itvnq 



(12) 



It follows from the properties of Q(p,uj) that the differ- 
ential rate (fT2l satisfies the detailed balance condition 



o-g(-w) = e " /T cr 9 (w). 



(13) 



Eq. (1131) implies that while cr q (u)) ~ cr 9 (— a;) at small 
energy transfers <C T, the differential rate is expo- 
nentially small at large negative ui. 

To gain further understanding of the differential rate 
at small momentum transfer, we consider first the regime 
of relatively high temperatures T » T s , when the inter- 
action in Eq. ([2]) can be neglected (except for very tiny 
energy transfers, see below). The correlation function in 
Eq. (fT2|) is then easily evaluated resulting in 



2 f 4 

= v^r- / II dk * fkjk 2 (fk 3 + i)(/ fe4 + 1) (14) 

- h) 6(£ki + 6c 2 - 6c 3 - 6c 4 + w )> 



i=l 

x <5(fci + fc 2 - fa 



where is the Bose distribution. (Eq. (| 14[) can also be 
derived by using Fermi's Golden Rule.) 

At T s -C T <C To the chemical potential is given by 



P = ~fJ>o, Po = T 2 /T Q < T. 



(15) 



At \u\ <C T the differential rate is dominated by pro- 
cesses in which both the initial and the final states of the 
two low-energy particles involved in a collision belong to 
the part of the spectrum with high occupation numbers: 
fkt « hi + 1 « r/(£ fe< +/i ) > 1. Evaluation of Eq. (HJ 
with this approximation results in 

a q (uj) = a 2 (T /^ 2 (To/TfF(\oj\/po), \u>\ « T. 

(16) 

The analytical expression for the function F(z) is some- 
what cumbersome. It is a monotonic function normal- 
ized as J Q F(z)dz — 1/8, with a power-law behavior at 



3 



z 1, F(z) = (2y/2/w)z 5 / 2 , and a logarithmic asymp- 
tote F(z) = (5/16tt 2 ) ln(8e- 7 / 5 /z) at z -> 0. 

The logarithmic divergence at cj — > in Eq. ([TB"]) comes 
from fci w fc2 ~ fc3 ~ ^4 in the integral over momenta 
in Eq. ([i"4"|) , and is an artifact of the free-boson approx- 
imation. The probability of scattering two bosons with 
close momenta k\ ~ k 2 in the initial state is suppressed 
at |&i — A&l ^ mc - (There is a similar suppression for the 
final states £3,4.) The logarithmic divergence in cr q (uj) is 
thus regularized at \u\ < uj 2 /T for T > T s . 

At T s <C T <C Tq and co 3> fJ.o, the main contribution 
to the integral in Eq. (fT4)) comes from | fex , 2 1 ^ V^Mo and 
|fc 3i 4| ~ yjrwjj >• I fex,2 j - Neglecting and £k lj2 in the 
arguments of the delta- functions in Eq. (|14p , we hnd 



(T q {uj) 



a 2 T 



<J2 



2vry^H (1-e— / 2T )' 



, M>Mo, (17) 



where 52 = 2 as appropriate for T 3> T s . Eq. (fl7|) ex- 
trapolates between Eqs. (fTTj) and (fl~6l) . 

With lowering the temperature, Eqs. (fl^j) - ([T5]) become 
inadequate when fj,Q (T) is of the order of the interaction 
energy per particle u> s , i.e., at T ~ T s . At T -C T s , 
however, the Bogoliubov approximation for the local den- 
sity operator becomes applicable (l9j . In this approxima- 
tion, excitations of a ID Bose liquid are essentially free 
phonons (Bogoliubov quasiparticles) , described by the 
Hamiltonian H B = J2k £ kb\b k with e k = y/£,k(£,k + 4w s ). 
In terms of phonons, the density operator has the form 
p{x) = n + EkM n ^/ L£ k) 1/2 ( b k + bl k )e ikx , where L is 
the size of the system. Using this representation, evalu- 
ation of Eq. ([T2"T) is straightforward and yields 



a q (w) 



647r v /27^Vl- e -"/ 2T 



u/uj s 



x [1 + (c/4 Ws ) 2 ]- 1/2 {l + [1 + ( W /4 Ws ) 2 ] 1/2 } 



(18) 

-3/2 



At \uj\ > oj 8 , Eq. ([18]) reduces to Eq. ([17]) with g 2 = 1 
appropriate for T <C T s , In fact, Eq. (TT71) is valid at any 
T <C To, provided that the energy transfer falls within 
the range max{/xo,a; s } <C <C £ g . For positive uj in 
this range, the validity of Eq. (fTT)) is due to the fact that 
the interaction has negligible effect [lj| on the final states 
of the colliding particles (£fe 3 « £fc 4 « ui/2 ^> lu s ). The 
applicability of Eq. (|17l) for negative w in the above range 
then follows from Eq. ([TB")) . 

We show the typical plots of the differential relaxation 
rate <j q {u) at T > T s and T < T s in Fig. [Q 

We turn now to the evaluation of the total relaxation 
rate, Eq. ([7]) . There are two contributions to the integral 
over oj in ([7]): 

r^r^ + fv (19) 

The first contribution, T^, comes from the high-energy 
"tail" of &q(uj), see Eqs. (fit))) and (TTT1) . This contribution 



T»T S 




FIG. 1: The differential inelastic scattering rate <T q {ui) at 
different temperatures (only |w| <C £, q domain is shown). The 
two plots correspond to Eqs. (|14[) an d HE}, evaluated at 
T = 4T S and T = T s /4, respectively, with T s = T /16. 



is independent of q and is given by 18 1 

a 2 T g 2 



r — 
1 \ — 



3^ 



(20) 



Note that, unlike in higher dimensions, Too depends on 
temperature via g 2 (T), see the discussion above. 



The second contribution in Eq. (1T91) . T q cx £ g 



1/2 



comes from the processes with a small energy transfer 
< max{T, oj s }. Using Eq. (fTBj) . we find 

f .=^(ff(if^ — »• 

(21) 

At lower temperatures we obtain, with the help of 
Eq. (HID, 



~ _ o^Tq fT s \ 1/2 fT \ 3/2 fT 
q ~ 16 



T. 



T 2 

-f < T « T s 



(22) 

Comparison with Eq. (|20|) shows that for not too large 
energies, £ g <C T s (Tq/T s ) 3 , the small momentum transfer 
contribution T q dominates the relaxation rate (|19[) in a 
broad temperature interval 



• t \V4/T \ 3 / 4 /T \V2/T \ 3 / 2 

Ht) (I) (I) '< 23 > 

which includes T = T S . At some temperature T m ax^ T s 
within this interval, the relaxation rate reaches its peak 
value r max = r g (T max ), see Fig. [2] By extrapolating 
the asymptotes ([21]) and (I22|) to the region T ~ T s and 
finding their intersection, we estimate T max « 1.6 T s , and 



3/4 



r„ 



0.16 a 2 ^^/^) 1 / 2 ^/^) 3 / 2 . (24) 



The actual values of T max and r max may differ from 
the above estimates only by numerical factors; finding 
these values would require a systematic description of 
the crossover regime T ~ T s . 

We now discuss briefly the feasibility of observing re- 
laxation by inelastic collisions in a system of cold atoms 
confined in a cylindrical trap. In this case the effective 



4 




FIG. 2: Sketch of the temperature dependence of the total 
inelastic relaxation rate. The dependence is nonmonotonic, 
with a maximum at T max ~ T s The dashed lines indicate the 
high- and the low-temperature asymptotes Eqs. (fl9"jl - (f22"]l . 



Hamiltonian ([I])-© can be derived explicitly, by projec- 
tion onto the lowest subband of transverse quantization. 
For a model in which the interaction in 3D is described 
by a pseudopotential V3r>(r) = 4n(a/m)S(r), where a is 
the s-wave scattering length [3j , and with the amplitude 
of radial zero-point motion a r = (ma; r )~ 1/ ' 2 ^> a (here 
uj r is the trap frequency), one finds foL [20j 



7 = 2a/na 2 , 



181n(4/3)(a/a r ) 2 



(25) 



The main limitation arises due to 3-body recombina- 
tion processes [2lj, absent in our model. The correspond- 
ing rate is T^ = f3n 2 g 3 /af 



2l|, where g 3 = (: p 3 :)/r 



Using Eqs. (J20J and ([25]), we find [18| 

Too/rfl = V92/93, V = 10.3a 4 / (m/3) 



(26) 



For 87 Rb (a = 5.3 nm, /3 = 3x 1CT 31 cm 6 /s [22|), we have 
r\ « 20. For a weak to a moderately strong interaction, 
7 < 1, the ratio 32/53 in Eq. f|26f) is of the order of 1 at 
all T, and T^/Tr ~ 10. 

In experiments with periodically colliding clouds of 
cold gases [llj, the 3-body recombination occurs all the 
time, while the scattering between the clouds takes place 
only during the collision itself (about one tenth of a pe- 
riod in [llj]). Therefore, the probability that during a pe- 
riod a particle participates in an inelastic collision event 
with a large energy transfer, and the probability that it 
participates in a 3-body recombination process are of the 
same order. Accordingly, inelastic scattering with a large 
energy transfer is difficult to detect unambiguously. 

Relaxation by the inelastic scattering with a small en- 
ergy transfer is effective when the interaction is weak, 
7 < 1 (indeed, the interval T*/T < T < T , see Eqs. 
(|2T|) and (l22t , disappears for large 7) . For the peak value 
of T q [see Eq. ([24])], we find 



(27) 



which for a fixed ratio £, q /T s diverges in the limit 7 — > 0. 

The maximum of T q is reached at T ~ T s . The con- 
dition of the observability of the inelastic relaxation, 
Tmax ^ and the condition for the high-energy 



quasiparticle to be outside the quasicondensate yet well 
within the lowest subband of transverse quantization, 
T s <C i q <C w r , can be satisfied simultaneously. For 
example, for 87 Rb the trap frequency w r /27r = 15 kHz 
and concentration n — 7/im~ 1 correspond to 7 = 0.2 
and T s = 120 nK. For ^ q /T s = uj r /£ q = 2.4, Eq. ([27]) 
then yields T max /Tn ~ 100. The parameters above are 



realistic with today's experimental technology [111, [2 



In conclusion, we evaluated the quasiparticle relax- 
ation rate in a weakly-interacting ID Bose liquid. Unlike 
in 3D, the rate is strongly momentum and temperature 
dependent, with a maximum at T ~ ns, where s is the 
sound velocity at T — 0. Our predictions can be verified 
in experiments with colliding clouds of cold atoms. 
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